Dinitrosyl formation as an intermediate stage of the reduction of NO in the presence 

of M0O3 
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We present first-principles calculations in the framework of density-functional theory and the 
pseudopotential approach, aiming to model the intermediate stages of the reduction of NO in the 
presence of MoO3(010). In particular, we study the formation of dinitrosyl, which proves to be an 
important intermediate stage in the catalytic reduction. We find that the replacement of an oxygen 
of M0O3 by NO is energetically favorable, and that the system lowers further its energy by the 
formation of (NO) 2. Moreover, the geometry and charge distribution for the adsorbed dinitrosyl 
indicates a metal-oxide mediated coupling between the two nitrogen and the two oxygen atoms. We 
discuss the mechanisms for the dinitrosyl formation and the role of the oxide in the reaction. 

[To be published in the Journal of Chemical Physics] 
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Reduction of nitric oxide (NO) has long been the focus 
of detailed studies, because of the role of NO as pollutant 
in the atmosphere. Three-way catalysts, based on Rh, Pt 
or PdO are used in automobile exhausts for this purpose; 
the effectiveness of such catalysts was found to increase 
when molybdenum oxides are addecU^. Recently, it was 
proposed that the NO reduction products N2 and O2 are 
formed nondissociatevely from NO via an adsorbed dini- 
trosyl species, which facilitates N-N coupling 3 ., since low 
temperature NO coupling proceeds through a dinitrosyl 
intermediate^. The geometry of the adsorbed dinitrosyl 
and the mechanism of N-N coupling remain unknown. 
To address this important question, we perform first- 
principles calculations based on density-functional theory 
and the pseudopotential approach, for various configu- 
rations involving nitrosyls and the MoO3(010) surface. 
The paper is organized as follows: First, we review the 
structure of M0O3 and describe the method we use for 
the calculations. Next, we present our results for the 
MoO3(010) surface with an O vacancy, a structure with 
an exposed Mo atom which can adsorb NO molecules. 
This is followed by detailed analysis of the configurations 
involving adsorbed NO and (NO)2- Finally, we discuss 
the dinitrosyl formation mechanisms and the role of the 
MoO3(010) surface as a catalyst for the reduction of NO. 
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FIG. 1: Structure of bulk Mo03. Mo atoms are shown as 
larger gray spheres, O atoms as white, smaller spheres. The 
three different types of O atoms (terminal, asymmetric bridg- 
ing, symmetric bridging) are indicated. 
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I. MOLYBDENUM TRIOXIDE 

Molybdenum trioxide is a layered material, its lay- 
ers being weakly bonded via van der Waals interac- 
tions. The space group is Qjf (bum), and the lat- 
tice is orthorhombic 5 with parameters 6 3.962, 13.858, 
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and 3.697 A. Each unit cell contains four M0O3 units. 
The lattice coordinates of the atoms are ±(x,y, 4) and 
±(| - x, \ + y, i), with (x, y) equal to (0.086, 0.099) for 
the Mo atoms , and (0.086, 0.250), (0.586, 0.099), (0.586, 
0.431) for the three O atoms surrounding each Mo atom, 
respectively^. 

The M0O3 crystal is shown in Fig. 1. Macroscopically, 
the material forms "small, thin, lustrous plates, parallel 
to (010)"^. This is revealed in the atomistic structure: 
It consists of bilayers parallel to the (010) plane, which 
are bonded through weak electrostatic interactions, with 
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the (010) surface being the easy cleavage plane of the 
crystal. Each bilayer consists of two sublayers of peri- 
odically arranged distorted MoOg octahedra. There are 
three structurally different O atoms. The asymmetric 
bridging oxygen is collinear with two Mo atoms and forms 
one long and one short bond with them. The symmetric 
bridging oxygen is located between the two sublayers of 
the bilayer and bonds to two Mo atoms of one sublayer 
with equal bond lengths and to one Mo atom of the other 
sublayer with an elongated bond. Finally, the univalent 
terminal oxygen is connected to one Mo atom forming 
the shortest Mo-0 bond in the system. The Mo-0 bond 
lengths in M0O3 as obtained from several experimental 
and theoretical studies are shown in Table [i] 



II. THE CALCULATION 

All calculations reported in this work were performed 
via the High-performance-Fortran Adaptive Real-space 
Electronic Structure (HARES) package 10 . We use 
Density-Functional Theory (DFT) in the Local Spin Den- 
sity Approximation (LSDA) 11,12 .. The Kohn-Sham va- 
lence electron wavefunctions are represented in a real- 
space orthogonal grid with a spacing of 0.13 A, which 
is equivalent to a cut-off energy of 150 Ry in a plane- 
wave basis calculation The ionic cores and their inter- 
action with valence electrons were taken into account 
through the soft Trouper-Martins^ pseudopotentials, in 
the separable Kleinman-Bylander formic. The Laplacian 
is expanded to second order terms. For the exchange- 
correlation functional we use the results of Ceperley and 
Alder as parametrized by Perdew and Zunger— . The 
electronic density was updated during the self-consistent 
loop using a modified Broyden mixer—. 

We model the surface of M0O3 by a single bilayer per- 
pendicular to the (010) direction. This approximation 
is not expected to affect the results, due to the weak 
inter-layer coupling. The calculated Mo-0 bond lengths 
for the single-bilayer system are shown in Table [I] As 
experimental data have discrepancies, we compare our 
results to the average of the two sets (last column of 
Table [|J. The agreement to experiment is within 2%. 
Our results are also in agreement to the calculation of 
Yin et al&. Moreover, vibrational frequencies calcu- 
lated within the single bilayer approximation^ are also 
in very good agreement with experiment. This justifies 
our choice to use a simplified unit cell, while establish- 
ing the adequacy of DFT and HARES to correctly de- 
scribe bonding in the material under study. The sur- 
face unit cell we chose has (2 x 2) periodicity relative to 
the ideal structure. This is sufficient to simulate isolated 
O vacancies or adsorbed NO, due to the local charac- 
ter of the substrate-adsorbate interaction and the lack 
of reconstruction effects for this surface. The Brillouin 
Zone was sampled by the T point only, which is adequate 
due to the large size of the unit cell. The atomic de- 
grees of freedom were relaxed with the Broyden-Fletcher- 



Goldfarb-Shanno (BFGS) methodic until the calcu- 
lated Hellman-Feynman forces were smaller in magni- 
tude than 0.005 a.u. We used the symmetry of the 
problem to eliminate some electronic as well as ionic de- 
grees of freedom where appropriate (see below). We have 
checked these computational parameters against previous 
calculationsiSii^ of similar systems to make sure they give 
well-converged results. 



III. TERMINAL O VACANCY AND NO 
ADSORPTION 

An O vacancy in M0O3 is formed when a terminal 
oxygen atom is missing. As discussed in Section [I] the 
terminal O is univalent and forms the shortest Mo-0 
bond in the bulk; this bond is also the strongest, with 
a calculated^ bond order of 1.93. Accordingly, it is ex- 
pected that the formation energy of the vacancy will be 
high. Indeed, our calculation yields a formation energy of 
2.95 eV/vacancy, assuming that the O atoms that leave 
the surface form O2 molecules (see Table llTlll . The (2 x 2) 
unit cell corresponds to a 0.25 ML coverage of vacancies, 
which is low enough to ensure that the interactions be- 
tween vacancies are negligible. The formation of the O 
vacancy is accompanied by small relaxations of atoms 
around the missing terminal O atom. The relaxed geom- 
etry is shown in Fig. 2(a). The main effect is an inward 
displacement of the exposed Mo atom, which results in 
shrinking the Mo-symmetric oxygen bonds from 1.96 A to 
1.91 A for O in the same sublayer and from 2.33 A to 2.11 
A for O in the adjacent sublayer. None of the atoms has 
any significant displacement parallel to the (010) plane; 
we thus find no evidence of surface reconstruction due to 
vacancy formation. 

The MoO3(010) surface with a terminal oxygen va- 
cancy, shown in Fig. 2(a), has many characteristics of a 
good catalytic platform: The exposed Mo atom promises 
to be chemically active and has the ability to bond to 
more than one adsorbate atoms, since it has lost a dou- 
ble bond. The remaining terminal O atoms around the 
vacancy site form a cage structure that can enhance 
coupling between the adsorbates. In ideal M0O3, the 
calculated^ population for a terminal O is -0.4 electrons, 
so the remaining terminal O atoms are expected to repel 
the adsorbates, creating an implicit attraction between 
them. 

The first step in the catalytic reduction of NO is the 
adsorption of a NO molecule. The binding configuration 
is shown in Fig. 2(b). NO is bonded to a Mo atom and 
substitutes a terminal O. The molecule is perpendicular 
to the surface in the lowest energy geometry, rendering 
the dissociation of nitrosyl on the surface implausible. 
The adsorption of NO on M0O3 with an O vacancy re- 
stores the positions of the surrounding atoms to their 
positions in the ideal MoO3(010) surface: the lengths of 
the Mo-N and N-0 bonds are 1.91 and 1.14 A, respec- 
tively. The N-0 bond has very similar length to the bond 
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TABLE I: Bond lengths (in A) between Mo and O in M0O3 as calculated in the present work compared to other calculations 
(a) by Chen et air- and (b) by Yin et al& and to experimental data obtained from (c) the atomic positions given by Wyckoff 6 
and and (d) the work of Kihlborg 9 as given in Refi— . The symbols Ot, O s and O a denote the terminal, symmetric bridging and 
asymmetric bridging O atoms respectively. 



Bond 


Bond Length 




This work 


Theory (a) 


Theory (b) 


Exp 


(c) 


Exp. (d) Exp. 


(avg) 


Mo-Ot 


1.77 


1.67 


1.76 




1.82 


1.67 


1.75 


Mo-O s 


1.96 


1.92 


2.02 




1.93 


1.95 


1.94 




2.33 


2.30 


2.28 




2.32 


2.33 


2.33 


Mo-Ot 


1.87 


1.76 


1.81 




1.89 


1.73 


1.81 




2.17 


2.19 


2.25 




2.08 


2.25 


2.17 



of a free NO molecule. We find this bond to be 1.15 A, 
in excellent agreement to experiment. 

The binding energy for NO to the MoO3(010) surface 
with an oxygen vacancy is 2.98 eV/molecule. The en- 
ergy difference between the MoO3(010) surface and the 
same surface with NO substituting a terminal O is —0.03 
eV/molecule, assuming that the O atoms leaving the sur- 
face form O2 molecules. The system lowers its total en- 
ergy by the substitution of terminal O by NO, so this 
configuration must be an important step in the NO re- 
duction. The mechanism of the substitution and the for- 
mation of O2 gas is not known, but it seems unlikely 
that the system passes through the O vacancy configura- 
tion, since the barrier for this path would be at least 2.95 
eV, a number which would forbid the process under ordi- 
nary conditions. The transition state for the substitution 
could involve intermediate formation of NO2 and/or O3. 

NO adsorption could also take place after a vacancy on 
the surface has been formed by an independent mecha- 
nism, for example after a CO molecule has hit the surface 
and CO2 has been formed. Another similar mechanism 
would be the formation of H2O from gas-phase H2 and a 
terminal O atom. The exact final product of the reduc- 
tion of MoO 3 (010) to MoO 3 [v](010) cannot be ad hoc de- 
termined, as it depends on the experimental conditions. 
However, the final state of the O atom that leaves the 
surface to create the vacancy would not affect the pic- 
ture, as the energetics would not change dramatically. 
We therefore choose to use O2 as a typical product of a 
reduction reaction. 



IV. DINITROSYL CHEMISORPTION 

To study the energetics of the dinitrosyl adsorption, we 
consider two characteristic configurations: In the first, 
the (NO)2 plane is parallel to (001) or the plane defined 
by the Mo atom and the symmetric bridging O atoms. In 
the second configuration, the (NO)2 plane is parallel to 
(100), or a plane perpendicular to the previous one which 
contains the Mo atom and the asymmetric bridging O 
atoms. The relaxed geometries for both configurations 



TABLE II: Bond lengths (in A) for the various structures 
discussed in the text. Ot, O s and O a stands for the terminal, 
symmetric bridging and asymmetric bridging O, respectively. 
The Mo-O s bonds for Mo and O s belonging at the same or 
adjacent sublayers, as well as the short and long Mo-O a bond 
are shown, (s) and (a) indexes for adsorbed (NO)2 refer to 
adsorption with the molecule parallel to the symmetric and 
asymmetric bridging O planes, respectively. 



Structure 


Mo - Ot 


Mo-O s 


Mo-Oa 


Mo-N 


N-O 






same 


adj. 


short 


long 






Ideal 


1.77 


1.96 


2.33 


1.87 


2.17 






O vacancy 


1.77 


1.91 


2.11 


1.90 


2.06 






ads. NO 


1.77 


1.95 


2.23 


1.86 


2.15 


1.91 


1.14 


ads. (NO) 2 (s) 


1.78 


1.98 


2.36 


1.89 


2.16 


2.00 


1.14 


ads. (NO) 2 (a) 


1.78 


1.96 


2.29 


1.90 


2.18 


1.95 


1.15 














2.35 


1.13 



are shown in Fig. 3. 

In the first reflection symmetry with respect to 

the (100) plane was imposed, following the experimental 
results obtained on oxidized Mo(110)^. Removal of the 
reflection symmetry does not result in further relaxation 
or significant lowering of the total energy, indicating that 
in the lowest energy geometry the two NO molecules of 
the dinitrosyl are identical. By contrast, (NO)2 adsorbed 
parallel to the (100) plane is found to be strongly asym- 
metric, with one of the two NO molecules closer to the 
Mo center than the other. This is a direct consequence 
of the asymmetry of the underlying oxygen atoms. 

The bond lengths for both configurations, together 
with the bond lengths of the previously discussed struc- 
tures are summarized in Table [H] For adsorbed dini- 
trosyl in the Mo-symmetric O plane, the Mo-N bond is 
2.00 A long, 5% longer than in that of a single adsorbed 
NO. The Mo-N bond is weaker for adsorbed (NO) 2 as 
the bonding electrons of Mo are shared by two N atoms. 
The N-0 bond length is 1.14 A, identical to that of ad- 
sorbed NO. In the second configuration, having dinitrosyl 
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(b) 

FIG. 2: Relaxed geometries for MoO3(010) (a) with a termi- 
nal vacancy and (b) with an adsorbed NO molecule. Grey, 
white and black spheres correspond to Mo, O and N atoms, 
respectively. 



in the Mo-asymmetric O plane, the Mo-N (N-O) bonds 
have lengths 1.95(1.15) and 2.35(1.13) A, for NO above 
the long (2.18 A) and short (1.90 A) Mo-asymmctric O 
bond, respectively Both N and O of the adsorbed ni- 
trosyls are more strongly bound above the less strongly 
bonded asymmetric O. 

The asymmetric configuration, with (NO)2 parallel to 
the Mo-asymmetric O plane, has a slightly lower energy 
than the other, by 0.10 eV per (NO) 2 molecule. The 
binding energy to MoO3(010) with an oxygen vacancy is 
3.80 eV and 3.90 eV per molecule for the symmetric and 




(b) 

FIG. 3: Relaxed geometries for dinitrosyl adsorption (a) in 
the symmetric O plane and (b) in the asymmetric O plane. 



asymmetric configuration, respectively. The energy gain 
when (NO)2 substitutes a terminal O of MoO3(010) with 
formation of 2 gas is 0.85 eV and 0.95 eV per molecule 
for the two configurations. 

The large lowering of the total energy by 0.90 eV on 
average per molecule for (NO)2 adsorption, compared to 
0.03 eV per NO for the NO adsorption, indicates that 
bonding between the two parts of the dinitrosyl has taken 
place. This bonding is evident in Fig. 4, which presents 
a contour plot of the total electronic density minus a su- 
perposition of atomic electronic densities on the dinitro- 
syl plane, which coincides with the symmetric bridging 
O plane. According to our choice of shading, electronic 
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FIG. 4: Contour plot of total electronic density minus a su- 
perposition of valence atomic charge densities in the Mo-N-O 
plane for adsorbed (NO)2 parallel to the symmetric O plane. 
The background shade corresponds to zero; darker shades cor- 
respond to negative values, lighter shades to positive values. 

charge has moved from the darker to the lighter regions, 
relative to the charge distribution of a superposition of 
atomic densities. The white T-shaped clouds inside the 
slab correspond to symmetric bridging O atoms, while 
those at the bottom to terminal O atoms. Small gray 
spots in the middle of black regions correspond to Mo 
atoms. The two bright vertical complexes correspond to 
the two NO molecules. As the electronegativities of the 
involved elements dictate, charge has moved from the Mo 
atoms to N and O atoms, with O atoms having slightly 



more charge concentrated around them compared to N 
atoms. 

We observe an apparent electron sharing between the 
two N atoms, manifesting the N-N coupling due to the 
presence of the catalytic surface. A similar, but weaker, 
coupling exists between the O atoms of the dinitrosyl. 
The absence of bonding charge between the N atoms of 
the adsorbate and the surrounding terminal O atoms im- 
plies a repulsion between them. It is this repulsion that 
leads to the attraction between the N atoms. 



V. DISCUSSION 

The reactions leading to dinitrosyl formation, starting 
either with gas phase NO or with adsorbed nitrosyl, are 
summarized in Table lTTTl Although the adsorption of two 
gas NO molecules to form a dinitrosyl is exothermic by 
0.90 eV/molecule (this number is the average between 
the two characteristic geometries of the dinitrosyl dis- 
cussed before), the process is expected to have a large 
energy barrier due to the repulsion of the NO molecules, 
which, having an electrical dipole moment, would prefer 
to have an anti-parallel configuration. The combination 
of two adsorbed nitrosyls is also rejected as a mecha- 
nism for dinitrosyl formation, as the vacancy left behind 
makes the process energetically costly. Indeed, as shown 
in Table IIIII this reaction is endothermic by more than 
2 eV/molecule. The alternative is to fill the O vacancy 
by an O from the environment. In this case, the dinitro- 
syl formation is exothermic by 0.81/0.91 eV/molecule. 
This process, being the opposite of the NO adsorption, 
would involve the same transition state and consequently 
requires overcoming almost the same energy barrier. 

A lower energy barrier, with almost the same energy 
gain, is possible when a gas-phase NO molecule binds to 
an already adsorbed one (last reaction of Table ITTl*|) . The 
energy barrier for this case has to be lower than in the 
previous one, since this process involves no breaking of 
bonds, while the terminal oxygens of M0O3 surrounding 
the adsorbed NO will attract the N atom of the gas- 
phase NO. The calculated lowest energy, together with 
the expected low barrier, support the presumption that 
this is the dominant dinitrosyl formation mechanism. 

Combining the previous results, the proposed mecha- 
nism for the reduction of NO is described in the following 
series of reactions: 

2NO(g) -» NO(g) + NO(a) + ±0 2 + 0.03 eV 
-> (NO) 2 (a) + ±0 2 + 0.87 eV 
-> N 2 + 2 + 1.00 eV. 

A gas NO molecule is exchanged with a MoO3(010) 
terminal O, and the system lowers its energy by 0.03 
eV. Alternatively, as mentioned in Section IlIII a vacancy 
could be created by the reaction of the surface with some 
other gas-phase molecule; in that case the products of the 
reaction would not be necessarily include 2 . Another 
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TABLE III: Energetics of the reactions discussed in the text, as calculated in the present work. MoC>3(010)[v] refers to a 
MoO 3 (010) surface with a terminal O vacancy; (g) stands for gas-phase molecule and (a) for adsorbed molecule. The energies 
for the reactions involving adsorbed dinitrosyl are the averages of the two considered configurations mentioned in Section lT¥l 

Vacancy formation: 



MoO 3 (010) -> MoO 3 (010)[v] + §0 2 (g) -2.95 eV 
NO adsorption: 

MoO 3 (010)[v] + NO(g) -f [MoO 3 (010)[v] + NO(a)] + 2.98 eV 
MoO 3 (010) + NO(g) -> [MoO 3 (010)[v] + NO(a)] + |0 2 (g) +0.03 eV 
(NO) 2 adsorption: 
MoO 3 (010)[v] + 2NO(g) -> [MoO 3 (010)[v] + (NO) 2 (a)] + 3.85 eV 
MoO 3 (010) + 2NO(g) -> [MoO 3 (010)[v] + (NO) 2 (a)] + |O a + 0.90 eV 
(NO) 2 formation from adsorbed NO: 
2[MoO 3 (010)[v] + NO(a)] [MoO 3 (010)[v] + (NO) 2 (a)] + MoO 3 (010)[v] -2.09 eV 
2[MoO 3 (010)[v] + NO (a)] + |0 2 (g) -> [MoO 3 (010)[v] + (NO) 2 (a)] + MoO 3 (010) +0.86 eV 
[MoO 3 (010)[v] + NO(a)] + NO(g) [MoO 3 (010)[v] + (NO) 2 (a)] + 0.87 eV 



NO molecule binds to the already adsorbed one, lowering 
the energy of the system by about another 0.84 eV. The 
next step is either desorption of the two O atoms first (as 
O2) followed by desorption of the two N atoms (as N2), 
or desorption of one O atom first and then of the N2O 
molecule. In both cases, an O atom fills the O vacancy 
that is left behind so that the catalyst is left unchanged at 
the end of the process. The final state of the system, with 
N2 and O2 molecules is^i 1.90 eV lower in energy than 
the initial state of two NO molecules. The contribution of 
the catalyst in the process is twofold: first, it can reduce 
the barrier for NO reduction, and second, it makes it 
much more probable for two NO atoms to come close to 



each other and react. 
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